Critical fluctuations in a soliton formation of attractive Bose-Einstein condensates 
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We employ mean-field, Bogoliubov and many-body theories to study critical fiuctuations in the 
position and momentum of a Bose-Einstein condensate whose translation symmetry is spontaneously 
broken due to attractive interactions. In a homogeneous system, the many-body ground state 
of the symmetry-preserving Ifamiltonian is very fragile against superposition of low-lying states, 
while the mean-field theory predicts a stable bright soliton which spontaneously breaks translation 
symmetry. We show that weak symmetry-breaking perturbations cause the translation-symmetric 
many-body ground state to cross over to a many-body bright soliton. We argue that the center- 
of-mass fluctuations in the soliton state arise primarily from the depletion of the condensate to 
translation modes. We develop an extended mean-field theory to analytically reproduce these results 
obtained by the exact diagonalization method. 

PACS numbers: 03.75.Lm, 03.75.Nt, 03.65.Ta 
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I. INTRODUCTION 

Ultracold matter waves offer the possibility of ob- 
serving quantum-mechanical fluctuations associated with 
low-dimensional many-body effects and phase diffu- 
sions 0,111111,13. In Bose-Einstein condensates (BECs) 
of dilute gases, various techniques of controlling exper- 
imental parameters enable us to study a rich variety of 
phase transitions. A standard way to study a ground- 
state phase of the condensate is to introduce an order 
parameter with spontaneously broken symmetry. On 
the other hand, quantum fluctuations play a crucial role 
in the emergence of the order parameter, and the re- 
lated issues have been widely studied in the systems of 
BECs mmillillllllllll- Furthermore, since 
atomic condensates are on the mesoscopic scale, we may 
expect interesting finite-size effects which are absent in 
the thermodynamic limit. 

The exact many-body ground state of a one- 
dimensional (ID) system with translation symmetry does 
not exhibit the off-diagonal long-range order in thermo- 
dynamic limit. The ground state has the translation 
symmetry of the original Hamiltonian, but is extremely 
fragile against localization to a bright soliton when the 
interaction is attractive. By investigating the many-body 
energy spectrum of this system, we have found that there 
exists a large number of low-lying quasidegenerate many- 
body eigenstates, and the localization to a soliton is due 
to the superposition of these states In the thermo- 
dynamic limit, the energy gap between the ground and 
low-lying states vanishes, and therefore the localization 
occurs no matter how small the symmetry-breaking per- 
turbation is. Such a localized soliton 15] can be well 
described by the Gross-Pitaevskii (GP) equation. How- 
ever, in a mesoscopic system, an energy gap between 
the ground and low-lying excited states remains finite, 
and therefore we can expect a crossover between the 
state with the translation symmetry and the soliton state 



as a function of the strength of interaction. In such a 
crossover regime, quantum fluctuations are expected to 
be large and the localization process is nontrivial. 



The aim of this paper is to examine the critical fluc- 
tuations in a soliton formation of an attractive BEG. We 
consider two cases for the breaking of translation symme- 
try, that is, symmetry-breaking potentials and quantum 
measurements. We tackle this problem using the numer- 
ical diagonalization of the Hamiltonian and an extended 
mean-field theory. In particular, we explicitly calculate 
the uncertainty relation between the center-of-mass po- 
sition and momentum, and show that the inclusion of 
center-of-mass fluctuations of the soliton is crucial for un- 
derstanding the nature of the many-body ground state. 
These fluctuations emerge as a depletion of the conden- 
sate, and decrease with increasing the magnitude of the 
symmetry-breaking potential or the number of the mea- 
surements. 



This paper is organized as follows. In Sec.^ we intro- 
duce a model of a ID attractive Bose system with the pe- 
riodic boundary condition, and briefly review the many- 
body energy spectrum and ground-state properties of the 
uniform system. In Sec. IIIII we introduce as symmetry- 
breaking perturbations a symmetry-breaking potential, a 
double-well potential, and quantum measurements, and 
discuss critical fluctuations in the position and momen- 
tum of symmetry-broken states. In Sec. IIVI we construct 
the Bogoliubov ground state that incorporates the ef- 
fect of center-of-mass fluctuations due to the presence of 
the infinitesimal potential. We also propose an extended 
mean-field theory to describe the effect of quantum mea- 
surements on the uncertain relation between the position 
and the momentum. In Sec.0 we summarize the main 
results of this paper. 
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FIG. 1: (Color online) Schematic illustration of a quasi-one- 
dimensional torus. 



Figure [21 shows the low-lying many-body spectrum of 
the Hamiltonian in which extensive rearrangement 
of the distribution of the eigenstates is seen to occur at 
gN ~ 1 due to the quantum phase transition. We see 
that in gN > 1 the spectrum is classified by two indices: 
that is, the band index a and the angular-momentum 
index C characterize each level as 



H\C),^E,,c\C)a, (3) 



where 



II. QUANTUM PHASE TRANSITION IN AN 
ATTRACTIVE BEC 



C). 



(4) 



A. Many-body eigenstates 

We consider a system of weakly interacting bosons on a 
quasi- ID torus with radius R as schematically illustrated 
in Fig.^ Throughout this paper the length, angular mo- 
mentum, and energy are measured in units of i?, h, and 
fi^/{2TTR), respectively. The Hamiltonian of the system 
is given by 



H 



de 



(1) 



where 9 is the azimuthal angle, ip is the boson field op- 
erator obeying the periodic boundary condition ipi^) = 
'0(0-|-27r), and the interaction is assumed to be attractive 
{g > 0). Although the Hamiltonian Q is exactly solv- 
able and the exact solution is analytically obtained by the 
Bethe ansatz |16). in principle, the solution is intractable 
for a large number of atoms. On the other hand, the con- 
figuration interaction method can evaluate the ground- 
state energy rather accurately fl^. Experimentally the 
quasi- ID ring-shaped BEC has recently been realized in 
a circular magnetic trap 0. 

We employ a numerical diagonalization method to in- 
vestigate many-body eigenstates of the Hamiltonian 
by restricting the Hilbert space to that spanned by the 
angular-momentum I = 0, ±1 states. It is possible to ex- 
tend the calculation to including the I = ±2, ±3 states. 
However, all results obtained for I = 0, ±1 remain qual- 
itatively unchanged, and we thus restrict ourselves to 
I = 0,±1 unless otherwise stated (see Sec. IHII3ll . The 
Fock-state bases are thus written as |n_i,no,ni), where 
ni is the number of atoms with angular momentum I. 
The number of atoms and the total angular momentum 
are given by J2i = -^j J2i ^"-^ — respectively. 
The field operator is thus expanded as 

1 



(2) 



where c; is the annihilation operator of a boson with an- 
gular momentum I. 



The energy levels E^- c in band a are distributed accord- 
ing to E„^o < E„^±i '< E„^±2, . For < giV < 1, on 
the other hand, some Eo-^'s are almost degenerate. To 
understand the nature of the rearrangement at gN ~ 1, 
we first study the eigenstates in the absence of interac- 
tion gN = 0, where the eigenstates can be described 
by Fock states rto, ni). We define the Jth state 

(J = 0, 1, 2, . . . , ) as the one in which J atoms are ex- 
cited. For examples, the J — state corresponds to the 
noninteracting ground state given by 



J = 0: |0),=o = |0,iV,0), 



(5) 




FIG. 2: (Color online) Many-body excitation spectrum 
Ea.c — Eofi of the Hamiltonian H with A'^ — 200, where £ 
is the angular-momentum index and a is the index of the 
bands that appears for gN > 1. We plot only low- lying states 
with indices cr = 0, 1, and 2. The state \C)cr denotes the 
many-body eigenstate whose total angular momentum is £, 
and band index is a. 
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FIG. 3: (Color online) Fano factor ^(6*) in Eq. of the 
many-body ground-state wave function as functions of the 
strength of the interaction gN and azimuthal angle 6 with 
iV = 200. 



and the excited states are given as 



J = 
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=0 = |3,iV- 


3,0), 



(6) 



The Jth state is hence (J-l-l)-fold degenerate. The index 
cr, which is the band index in 17 TV > 1, corresponds to the 
number of Z — ±1 pairs in the Fock state for each Jth 
state at gN = 0, i.e., cr = ( J - \C\)/2 = 0, 1, 2, ... , . 
For < gN < 1, the levels are characterized by J alone, 
since the degeneracy with respect to cr at gN = is 
maintained. 

The excitations in < gN < 1 have substantial en- 
ergy gaps until the critical point gN ~ 1 is reached. On 
the other hand, for gN > 1, the energy difference be- 
tween the ground and the first excited states scales as 
We thus expect that the ground state for gN > 1 
is vulnerable for large N, and perturbations of the order 
1 /N can cause drastic reconstruction of the ground state. 
The emergence of these quasidegeneracies in the energy 
spectrum may be regarded as a precursor of symmetry 
breaking, and interesting effects such as the enhancement 
of the condensate fraction can be expected as a con- 
sequence of symmetry breaking. 



B. Ground-state properties 

We study properties of the ground state \C — 0)cr=o 
as a function of the strength of interaction gN. Ac- 
cording to the standard definition of the BEG j^^, the 
condensate fraction and the condensate wave function 



for a many-body state \'^) are given in terms of the 
maximum eigenvalue Am and the corresponding eigen- 
function of the reduced single-particle density matrix 
p{0',e) = (^'|7/'t(6»')T^(6»)|^')/iV. If the maximum eigen- 
value Am is on the order of unity and the other ones are 
on the order of there exists a usual single conden- 
sate described by the corresponding eigenfunction. For 
the case of gN > 1, however, the ground state does not 
fall into this category of BEG. By calculating the eigen- 
value of p for the ground state as a function of gN, we 
find that Am is on the order of 1 for gN < 1, but for 
gN > 1 there appears more than one eigenvalue on the 
order of unity. The ground state is thus a conventional 
single condensate for gN < 1, and fragmented for 
9N>1. 

Since we consider the Hamiltonian with translation 
symmetry, the expectation value of the number density 
of particles n{9) = Np{9, 9) of the ground-state wave 
function is constant for < < 27r. In contrast, the 
two-body correlation function. 



g<^^He,e') 



(7) 



is found to deviate greatly from unity for gN > 1, while 
it is almost unity for gN < 1 [2l|. We show in Fig. |31the 
Fano factor of the ground state defined by 



m ^ 



{n^{e)) - {n{e)f 



n{e) 

= l + n{9) \g^'^\9,9) - 1 



(8) 



The large deviation of F{9) from unity for gN > 1 indi- 
cates that the number density of particles has quantum 
fiuctuations, which may also be regarded as a precursor 
of formation of the broken-symmetry state. 



III. SYMMETRY-BREAKING 
PERTURBATIONS 

In this section, we investigate on how the many- 
body ground state responds to infinitesimal symmetry- 
breaking potentials or quantum measurements. 



A. Infinitesimal potential 

1. Superposition of low-lying states 

We employ the exact diagonalization method to obtain 
the ground state of the Hamiltonian 



Ki=H + eiVi, 

r27r 



(9) 



where Vi = J^"" d9'iP^9) cos 9ip (9) is a symmetry- 
breaking potential. Since we discuss here only low- 
lying modes (cr = 0, £ ^ 0) close to the ground state 



(cr = 0, £ = 0), we simplify the notations of the eigen- 
states by omitting the index a. 

We consider the energy change in the many-body 
ground state as a function of the magnitude ei of the 
potential. Because of the degeneracy Ec — E^c with 
respect to the angular momentum, we assume that the 
broken-symmetry ground state 1^'^'^^'') of the Hamilto- 
nian Ki is described by 



I* 



(^i)\ 



-iLe 



/3o|0)-t-^/3z:(|/:> 



£>0 



, (10) 



where the coefficients (3c satisfy the normalization con- 
dition, J2c\M^^l- 

In Fig. 21 we show the energy change in the ground 
states 



AE„ 



^ {¥;^^\k,\¥;^^)-{o\m (<o), (ii) 



multiplied by N, as a function of 
Kl EE eiN'^. 



(12) 



The slope of the ki dependence of AEgs changes near 
Ki ~ 1. The Ki dependence of AEg^ for < ki < 1 
can be estimated by perturbation theory. The first-order 
correction to the ground-state energy is zero, and the 
second-order correction is given by 



£#0 



Eo-Ec 



(13) 



where the denominator and the numerator are on the 
order of A'^"^ and (eiN)^, respectively. The energy dif- 
ference AE'gs therefore depends on —efN^ = —k\/N. 
At the transition point ki ~ 1, the shift in the ground- 
state energy due to the symmetry-breaking perturbation 
is A|£'gs| ~ 1/-/V, which is on the same order of magnitude 
as the energy gap. As ki exceeds unity, the energy scale 
of the potential exceeds the energy gap E±i — Eq and the 
matter wave begins to localize by superposing the ground 
and low-lying excited states as shown in Fig. |S{a). The 
perturbation theory, therefore breaks down for ki > 1. 
The Ki dependence in this regime can be well described 
by the mean-field theory as shown in Sec. IIVB II In the 
absence of the symmetry-breaking potential, the forma- 



tion of the broken-symmetry state 
by an amount of 



costs energy 



mm 

= Y,{Ec-E^)\Pc? (>0), (14) 
c 

which gives an energy increase associated with the super- 
position of the low-lying state. 

The symmetry breaking must be associated with a 
Nambu-Goldstone mode, which plays the role of restor- 
ing the symmetry. In the present context, the Nambu- 
Goldstone mode should be the translation zero mode 
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FIG. 4; The amount of decrease in the ground-state energy 
given by Eq. pi|l in the presence of a symmetry-breaking 
potential for gN — 1.5, where = e\N^ . The superim- 
posed dashed line shows the result of the mean-field calcula- 
tion given in Eq. I|87|l . 



|$i^'^'') of the localized state l^'i^'^'') which is defined as 



Y,CPc{\C)-\-C)). (15) 



£>0 



The broken-symmetry state with the excitation to the 
translation mode is thus given by 



I* 



1 



(16) 



where C, denotes the relative amplitude of the excitation, 
and the energy difference associated with the excitation 
of the translation mode is given by 



1 



ci>^)|Xi|ci>rO), (17) 



where A^tr = \J {^t^^^t^^) ■ We note that Etr becomes 
zero in the limit of C — > 0. 

2. Quantum fluctuations in position and momentum 

Let us investigate the quantum fluctuations in the 
ground state l^'g^^^) as a function of the magnitude of the 
symmetry-breaking potential. The angular-momentum 
distribution is given by \l3c\'^ = \{C\'^'f^^)\^, which be- 
comes Gaussian-like for ki ^ 1 as shown in Fig. IH^a). 
In Fig. EJb) we plot the variance (A£)^ of the angular- 
momentum distribution with open circles. The devia- 
tion in Ki dependence of (A>C)^ in < ki < 1 arises 



because in this region the energy scale of the symmetrj 
breaking potential is smaller than the energy gap betwee 
the ground and the first excited states, and the angula: 
momentum fluctuation is suppressed. Once the matte 
wave is localized by the superposition of the form (|f 0| 
the fluctuation A>C of the angular-momentum distribr 
tion is proportional to kJ^^. 

It follows from the uncertainty relation A£A^?c.in. — 
that the center-of-mass fluctuation of the localized stat 



I* 



(ei) 



for Ki > 1 is given by 



-1/4 



(1? 



which is significant for small perturbations j22 |. I 
we will show that this ki dependence is we 



Sec UVBl 

described by the Bogoliubov ground state that takes int 
account this position fluctuation by the depletion of th 
condensate. 



3. Condensate fraction 



By using the many-body ground state I'^f^'') in th 
presence of the symmetry-breaking potential, and d: 
agonalizing the single-particle density matrix p(l, V) - 
(^'fl^^^|c|,c;|^'fl^^''), we obtain the three eigenvalues in th 
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FIG. 5: (Color online) (a) Density profile n(Q) of the many- 
body wave function and (b) the peak density as a function of 
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FIG. 6: (Color online) (a) Distribution of the angu- 

lar momentum £ in the many-body ground state l^*^^^') for 
gN = 1.5. (b) Width of angular-momentum distribution in 
the many-body ground state 1^''^^') (open circles) and that 
in the Bogoliubov ground state 1^''-^-') (solid line) given by 
Eq. ra . 



truncated bases. We define the number of depleted atoms 
N' of the many-body ground state as N' /N = 1 — A^^i-', 
where A*-'^^'* is the maximum eigenvalue of p. We plot 
N' /N in Fig.[71as a function of ki. For ki = 0, the max- 
imum depletion is obtained due to the fragmentation of 
the many-body ground state, and it remains constant for 
Ki < 1. The depletion N' suddenly begins to decrease at 
Ki ~ 1, and the ki dependence for ki > 1 is found to be 



-1/2 



The condensate fraction thus increases when the 



ground state begins to locahze. 



B. Case of a double- well potential 

We next discuss what happens when there exists a sym- 
metric double- well potential of the form V2{0) = cos 26*, 
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FIG. 7: Depletions of the condensate obtained by the di- 
agonahzation of the Hamiltonian (open circles), and by the 
Bogoliubov theory given in Eq. H57|l (solid line). 



and discuss the ground state of the Hamiltonian 



where 



1/2 



K2 = H + e2V2, 



dei^He)V2(9),jj{e). 



(19) 



(20) 



In the double-weU potential, a rich variety of experi- 
ments, e.g., diagnosing the correlation functions from 
interference and uncertainty relations, become possi- 
ble 

For an attractive condensate in a sufRciently deep sym- 
metric double- well potential, the ground state may form a 
Schrodinger's cat state, where the localized macroscopic 
states on the left and the right wells arc in a superposition 
state miillii, |23|. We diagonalize the Hamiltonian 
K2 including angular-momentum states I = 0, ±1, ±2 be- 
cause the symmetry of the potential requires the inclu- 
sion of the I = ±2 states. The solid curves in Fig.|Hlshow 
eigenvalues of the reduced single-particle density matrix. 
The presence of more than one large eigenvalue is a signa- 
ture of the Schrodinger's cat state. While the maximum 
eigenvalue Ai increases with increasing ei, it decreases 
slowly with increasing £2, and simultaneously the second 
maximum eigenvalue A2 begins to grow. These eigenval- 
ues Ai_2 eventually approach 1/2, corresponding to the 
states localized on the one or the other well. 

In the case of the double- well potential, the many-body 
ground state is also described by the superposition in the 
form of Eq. (|10l) . Figure |51 shows that the K2 dependence 

of the energy difference (^-^^'"^ |i^2|*i,^'^) - {0\H\0) is sim- 
ilar to that of the single-potential case shown as a dashed 
curve. This indicates that the superposition begins near 
k;2 ~ 1, since the energy at that point is of the order 1/N. 
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FIG. 8: (Color online) Eigenvalues of the reduced single- 
particle density matrix for the double-well potential (solid 
curves) as a function of K2 = £2N^ ■ Dashed curves show 
the results obtained for a single-well potential as a function 
of Ki = eiN'^. In the latter case only one large eigenvalue 
survives for ki > 1. 




FIG. 9: Ground-state energy as a function of K2 (solid curve) 
in the presence of a double-well potential measured from 
-Bo — (0|-ff|0). Dashed curve shows the results of a single- 
well potential. 



C. Quantum measurement 

We next consider quantum fluctuations caused by re- 
peated measurements. The issue of interference of the 
two independent BECs was discussed from the viewpoint 
of quantum measurements in Refs. [28l l29l |. Below we 
discuss an analogous process in which a fragmented con- 
densate makes a transition to a single condensate via a 
quantum measurement |l4l l30| . 



Quantum fluctuations in position and momentum 



The action of the quantum measurement at the posi- 
tion 9j relates the postmeasurement state j^'^-')) to the 
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FIG. 10: (Color online) (a) Cumulative number of counts 
for the center-of-mass position obtained from 2000 runs of in- 
dependent simulation of quantum measurements for gNi-ait ~ 
1.6 and Mnit ~ 300. (b) Center-of-mass fluctuations obtained 
by the numerical simulations (open circles) . The solid straight 



line obeys {A9c 



oc j and is drawn as a guide to the eye. 



FIG. 11: (Color online) (a) Ensemble average of the angular- 
momentum distribution \Pcf obtained from 2000 runs of the 
numerical simulation, and (b) its width (open circles) for 
gNinit = 1.6 and Mnit ~ 300. The solid line is given by 
Eq. ira . 



premeasurement state as 



^(^?,.)|*0-i)) 



;*o-i)|7/>t(0j)^(0j-)l^'O'- 



(21) 



where j is the number of measurements, and the initial 
state |\I>(J='^)) is taken to be the many-body ground state. 
The subsequent measurement position Oj^i after the jth 
measurement is probabilistically determined according to 
the density distribution nj{9) = Njp'^^^O.O) of the pre- 
measurement state, where Nj = A^init — j is the number 
of atoms. 

We perform 2000 runs of the numerical simulation of 50 
quantum measurements (|21|) {0 < j < 50) starting from 
the ground state with A^init atoms. By the measurements, 
the density distribution nj{0) is found to localize, and 



we denote the center-of-mass position of rij (6) by 0^,}^, . 

We find that as j increases, Oh m. converges to a certain 
position, which is random from run to run. We there- 

fore rotate the system so that 6*0.111. is located at the 
origin, and plot the distribution of 9c''.li. for 2000 runs of 
measurements in Fig. llOf a). Hence, this distribution rep- 
resents the center-of-mass fluctuation after the jth mea- 
surement. The variance of the distribution is plotted in 
Fig. llOf b') , which shows that 



(A0)2 cx 



(22) 



Because of the translation symmetry of the initial state 
the center-of-mass fluctuation is maximal before 
the measurement. The repeated measurement process 
reduces the position fluctuation, and eventually projects 
the translation-invariant state |vl'(j=")) onto a broken- 
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symmetry state whose center-of-mass localizes at a cer- 
tain position. 

Because the quantum measurement is a stochastic 
process, we consider the ensemble average of |/?£p = 
I ^£|v[/(j)^|2 over the 2000 runs of independent simulations. 
The result is shown in Fig- lllT a). We find from Fig. lllf b') 
that the distribution of the ensemble-averaged angular 
momentum obeys the Gaussian with its width given by 

(A/:)2 « j. (23) 

In the present ID system, the center-of-mass position 
^c.m. — X^fcLi ^k/N and the angular momentum C = 
obey the commutation relation 

%.n..,C]=i. (24) 

Their fluctuations A6'c.m. = ^^c.m. — (^c.m.) and A£ = 
C — (£) therefore obey 

A0c.m.A£ > i. (25) 

From Figs. [TUl and llll we found the uncertainty 

A6lc.m.A£ ~ 0.73, (26) 

which is larger than that of the minimum uncertainty 
state. 



2. Condensate fraction 

We study the change in the condensate fraction as a 
function of the number of measurements j. In Fig. 1121 we 
show the ensemble-averaged depleted fraction of the con- 
densate in the state for three different values of the 
initial number of atoms iVinit- The depletion decreases 
monotonically with increasing j . The change in the con- 
densate fraction with respect to the number of measure- 
ments j is not sensitive to the initial number of atoms 
iVjnit, and is determined by gMnit- Like a symmetry- 
breaking potential, the repeated quantum measurements 
change the fragmented condensate having a translation 
symmetry into a translation-symmetry broken single con- 
densate with reduced center-of-mass fluctuations. 



IV. BOGOLIUBOV AND EXTENDED 
MEAN-FIELD THEORIES 

In the previous sections we showed that center-of-mass 
fluctuations are significant when the symmetry-breaking 
perturbations are sufficiently small. We here analytically 
treat the many-body state with large center-of-mass fluc- 
tuations for both cases of the symmetry-breaking poten- 
tial and the quantum measurement. 




FIG. 12: Depletion of the condensate for gNh-iit = 1.6. Sym- 
bols (o,o,x) show the results of numerical simulations, and 
the solid curve shows a theoretical one in Eq. (Ifci9t . where j is 
the number of measurements. 



A. Mean-field ground state in a homogeneous ring 

By replacing the fleld operator with a c-number wave 
function as ip{9) (l){0), we obtain GP energy functional 

4° M = N r de{\S/cj,f - ngN\cj,f ). (27) 
Jo 

We employ a variational wave function 

(/)(0) = -^(ao+aie^'^+a-ie-'*'), (28) 
V 27r 

where ao,±i are complex variational parameters which 
obey the normalization condition |a_ip -I- |aoP -I- — 
1 and are to be determined so as to minimize Eq . For 
gN < 1, Eq is minimized when ao = l,Q;±i = 0, i.e., 
the condensate wave function is uniform, 0(0) — l/\/27r. 

fGP) 

For gN > 1, Eq is minimized when 

l 3gN + 2 , / 2(giV-l) ^^^^ 

where the global phase is chosen so that ao is real without 
loss of generality. The corresponding wave function 

(l>{0 - ^c.m.) = [ao + 2ai cos(0 - 0c.m.)] , (30) 

V 27r 

is a spontaneously broken-symmetry state. Note that all 
states described by Eq. H30(l are degenerate with respect 
to an arbitrary center-of-mass coordinate 9c.m.- The 
critical point gN = 1 of the quantum phase transition 
obtained by the present mean-fleld theory agrees with 
the exact analysis of the ID GP equation, and Eq. H29|) 
indeed corresponds to the dominant coefficients of the 
plane-wave expansion of the exact solution of the ID GP 
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equation for aN > 1 |2lLl31| . 



c.nO oc dn( ^^{9-e,.,r,.) 



oo 

1+E 



-ii{0-e^,m-) 



1=1 



cosh{TrlK'/K) 



-(31) 



where dn is one of the Jacobian eUiptic functions, 

K{m) = J^^^ dxj \/\ — m sin^ x is the complete elhptic 
integral of the first kind, and < m < 1. From the defi- 
nition of the GP mean-field theory, all atoms occupy the 
same single-particle state (/), and the condensate fraction 
is unity for all gN . For gN > 1, the ground state of the 
GP equation given in Eq. H3U|I is thus qualitatively dif- 
ferent from the many-body ground state |0)cr=o in both 
the condensate fraction and the translation symmetry. 



B. Infinitesimal potential 

1. Gross-Pitaevskii and Bogoliubov analyses 

We introduce the symmetry-breaking potential Vi{6) 
which is assumed to be infinitesimal, i.e., of the order 
of The GP energy functional in the presence of a 

symmetry-breaking potential Vi {9) = cos 9 is given by 



2ir 



Jo 

-7TgN\m\'] -(32) 

When we employ the same variational function of the 
form in Eq. (|28|) . e['^^'^ is minimized when the variational 
coefRcients ao and ai = a^i are given by 

2 2(g7V_i) 



S-gN 

7gN ' 7gNy'2{gN -l){3gN + 4) 



ei + 0{ei) 



(33) 
(34) 



where we assume that the center-of-mass position is lo- 
cated at 9c.m. = 0. The corresponding chemical potential 
and energy read 



-15gN 1 



14 



7gN 



%/2(8 - gN) 



^igN~l){3gN + 4) 



E- 



(GP) 



/N 



+ 



+ ^/2{gN ~l){3gN + 4) ei + 0(e^), (35) 
1 

28gN 
2V2 



15g^N^ + 16gN 



y/igN-l){SgN + A) + 0{el). (36) 



7gN 

From Eq. H3()|l we obtain 

7V(i?pP)-£;fP)) 

= -^^2{gN-l){3gN + A) + 0(e?),(37) 



which agrees with the ki dependence of NAEgs [see the 
dashed curve in Fig.^. 

In the case of uniform space (ei = 0), the translation 
mode can be excited without energy cost for gN > 1 21], 
which corresponds to the Nambu-Goldstone mode. In 
the presence of a nonuniform potential, the translation 
mode slightly costs energy. By solving the Bogoliubov- 
de Gennes equations 

-^^g-SlVl{9)-^I-2^rgN\^\^ -ngNcj)^ 
dl + eiVi{9)+ fi + 2TrgN\(l)\'^ TrgN(j)*^ 

(38) 

we obtain the eigenvector corresponding to the lowest- 
energy translation mode as 





u 




u 










V 




V 



u 




■ 1 + 




cx sin 9 




V 




1 - 



,3/2-) 

r3/2^| 



A = ±V^ 
where S is given by 



2(5, 



7ei 



/3.9iV-f 4 



, (39) 
(40) 

(41) 



2 + 5gN y 2gN - 2 
The corresponding quasiparticle operator is given by 

ho^uf + vf\ (42) 

where 



(43) 



The GP ground state whose center-of-mass position is 
located at ^c.m. is written as 



I* 



(GP)\ 



1 



ao4 + Q!i(c|e'^<^'"- + c^e"''^'^"-) |vac) 



N 



(44) 



and the derivative with respect to 6*0.111. gives 



d9, 



cx Tl* 



(GP) 



.=0/ 



(45) 



Thus, the operator T serves as the generator of the trans- 
lation of the GP ground state IvI'(C'P)). 



2. Quantum fluctuations in position and momentum 

We consider the Bogoliubov ground state. Since the 
translation costs little energy, the center-of-mass fluctua- 
tion is expected to be dominant in the Bogoliubov ground 
state. We therefore take into account the translation 
mode only. The Bogoliubov quasiparticle operator H42|l 
is written as 



6 = e- = ' 



f coshr + f^ sinhr,(46) 
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where the parameter r is related to the Bogohubov quasi- 
particle and hole amplitudes u and v by (cosh r, sinh r) oc 
(w, u). We evaluate the angular- momentum distribution 
in the Bogoliubov ground state 

|^(B)^ g§(f+^-f^)|^,(GP) \ 



.=0/ 



1 



/M 



:(aocJ, + \/2ai^^)^|vac),(47) 



where 



(48) 



Let us approximate the angular-momentum operator 

L^ni-n^i^S^f + Sf\ (49) 

as 

C = V2Nai{f + f''), (50) 

by the replacement of S with its mean-field value \/2N a\ 
under the assumption that the ratio of the condensate 
atoms to the excited atoms is sufficiently large. By using 
the relation (see Appendix^ for proof) 



e-5(*^'-*')|vac) 



1 



-e aU+'f I vac) 



coshr 

we obtain the Bogoliubov ground state as 



(GP) 



.=0/' 



(51) 



(52) 



where A/r is the normalization constant, and 
tanhr vju 



V 



1 — tanhr I — v/u 



5gN + 2 f2gN 



^ \ -1/2 -L 



14 



, +0(£|/^).(53) 

We thus find that the angular-momentum distribution in 
the Bogoliubov ground state is given by 



l(/:|*(B))|2 ^ 



1 



(54) 



v/27r(A/:(B))2 

where the width of this distribution is calculated to be 
(see Appendix IB|) 

2gN-2 



7gN 



with 



F 



7 



l3gN^ 



5gN + 2y 2gN -2 



(55) 



(56) 



The distribution of the angular momentum in j^'^^^) is 
in excellent agreement with the distribution coefficient 
|/3£|2 = \{C\^^g^^)\'^ in the many-body ground state ifTUI) 
which is obtained numerically by the diagonalization of 
the Hamiltonian K in Eq. 0. 



3. Condensate fraetion 

The condensate fraction of the Bogoliubov ground 
state is less than unity because of virtual particle-pair 
excitations mainly to the translation mode. The virtual 
excitations are the physical origin of quantum fluctua- 
tions in the position and momentum described above. 
The number of depleted atoms is calculated from Eq. H39|) 
as 



2Tr 



1 



(57) 



where 6 is given in Eq. (|41|l . The result is shown as the 
solid line in Fig. [3 which agrees well with the results of 
the exact diagonalization for ki > 1. 

We thus find that in the presence of the symmetry- 
breaking potential, the Bogoliubov ground state very well 
reproduces the energy, the angular-momentum distribu- 
tion, and the depletion of the condensate obtained by 
the exact diagonalization for ki >1. This is because the 
virtual excitations described by the Bogoliubov ground 
state lead to the center-of-mass fluctuations above the 
mean-field ground state. 



C. Quantum measurement 

1. Quantum fluctuations in position and momentum 

In this section the suppression of the center-of-mass 
fiuctuation via repeated quantum measurements shown 
in Fig. 1101 is investigated semiclassically by generalizing 
mean-field theory. We assume as an initial state super- 

(GP) 

position state of the GP solution ) with respect to 

the center-of-mass position ^c.m., 



(58) 



We assume that for j — 0, Ao{9c.m.) is a 
constant. The repeated measurements would se- 
lect a localized soliton at a certain 6'c.m.j '/'(^ — 
^'c.m.) = [q^o + 2ai cos(6' — 6c.m.)] /V2tt with probability 
^|(^c.m.)- If the first atom is detected at a position 6i, 
the postmeasurement distribution may be approximated 
by 



^l(0c.m.) 



1 



/2Tr 



[ao + 2aicos(0i-0c.m.)], (59) 



where ag and ai are given by Eqs. H29|) . Likewise, after 
the jth measurement, the distribution becomes 



1 

.) cx TT -= [ao + 2ai cos{9, - 0c.m.)] ■ (60) 
■ n v2'7r 



A. 



If we rotate the system so that Ocm. ^ for j — > oo 
(corresponding to the case of Fig. I10|) , the distribution of 
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6i should become [uq + 2ai cos 0]^/(27r). Equation 1601) 
is then simpUficd for j ^ 1 as 

j 

InAj (6*0.111.) = ^ In [ag + 2ai cos(9i — ^c.m.)] + const 

— j J d9n{9) In [ao + 2ai cos(6'i — 0c. m.)] + const 

~ -2jalel^^ + const + • • • . (61) 

The distribution function of the center-of-mass position 
^c.m. can thus be approximated as 



Aj{9c.m.) oc e 



(62) 



It foUows then that the center-of-mass position fluctuates 
according to 



A6'c.m. oc 



1 



(63) 



The angular-momentum distribution of the state (|58|l . 
can also be analytically calculated. Substituting Eq. 162|) 
in Eq. (|58|l . we find that the angular- momentum distri- 
bution function is given by 



|(/:|vi/(^'^)p(x-=4 



with the width 



(A£)2 = 2a?j. 



(64) 



(65) 



It follows from Eqs. (|63l) and (|65|l that during the re- 
peated measurements the fluctuations in the center of 
mass and the angular momentum are found to obey the 
uncertainty relation 



A0,.,i.A£ ~ 1. 



(66) 



Both Eqs. H63|l and (|65|l are in excellent agreement with 
the ensemble averages of the corresponding results ob- 
tained by the numerical simulation as shown in Figs. [TUI 
andlTTl 



2. Condensate fraction 

We showed in the previous section that the repeated 
measurements suppress the center-of-mass fluctuations. 
We show here that this leads to an enhancement of the 
condensate fraction. The reduced single-particle density 
matrix of l^'^f ) is obtained from Eqs. and (jHU 

as 



Pill') = {clci) 



at 



a 







-1//2 



aoaie 



-w/2 



, (67) 



where v = l/(8af j), and we assumed 



(68) 



for simplicity. The largest eigenvalue A^"* of p{l,l') is 
given by 



1 



+ ^(1 -al + a2e-2.)2_4c,2(i _ 2q,2)(i _ ^--y 



.(69) 



Figure El compares Eq. (|69|l with the results obtained by 
numerical diagonalization, where the quantitative agree- 
ment is found. Like a symmetry-breaking potential, the 
repeated quantum measurements also change the frag- 
mented condensate with a translation symmetry to a sin- 
gle condensate with a broken translation symmetry. In 
terms of the generalized mean-fleld theory developed in 
this section, the depletion of atoms is described by the 
amplitude Aj in the state (|58|l . the absolute square of 
which represents the semiclassical distribution function 
of the center of mass of the bright soliton. 



V. SUMMARY AND CONCLUSIONS 

We have investigated the critical fluctuations associ- 
ated with the formation of broken-symmetry state by 
explicitly introducing symmetry-breaking potentials or 
quantum measurements. 

In the absence of the symmetry-breaking potential, the 
many-body ground state, obtained by the diagonaliza- 
tion of the symmetry-preserving Hamiltonian, is found 
to be fragile against a formation of a localized state. The 
localized state corresponds to a superposition state of 
quasidegenerate low-lying modes. The localization is also 
reflected in the enhancement in the two-body correlation. 

In Sec. 1111 Al we have studied the many-body local- 
ized state in the presence of an infinitesimal symmetry- 
breaking potential. We have shown that the localization 
of the many-body ground state begins when the energy 
scale of the symmetry-breaking potential becomes of the 
same order of the energy gap between the ground and 
the first excited state of the symmetry-preserving Hamil- 
tonian, i.e., 1/N. We have evaluated center-of-mass fluc- 
tuations in the regime of the crossover between the state 
with the translation symmetry and the broken-symmetry 
state, and found that these fiuctuations are the origin of 
the depletion of the condensate fraction. In the presence 
of the double-well potential, we have shown that as the 
magnitude of the potential is increased, the ground state 
becomes a Schrodinger's cat state, i.e., the macroscopic 
superposition of two localized states. 

In a manner similar to the single-potential case, re- 
peated quantum measurements also cause the crossover 
from the uniform state to the localized state with sig- 
nificant enhancement in the condensate fraction. We 
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have shown that the center-of-mass and the angular- 
momentum fluctuations obey the uncertainty relation 
during quantum measurements. 
In the latter part of this paper, we have devel- 
oped an analytic method to treat the many-body states 
with quantum fluctuations. Starting from the broken- 
symmetry ground state of the GP equation, we have 
taken into account the translation mode of the soliton as 
a Bogoliubov fluctuation in the presence of an infinites- 
imal symmetry-breaking potential. We have found that 
the angular-momentum fluctuations and the condensate 
fraction of the Bogoliubov ground state very well agrees 
with those of the numerically obtained many-body local- 
ized state. 

The effects of the quantum measurement of the bright 
soliton on the localization have been studied by develop- 
ing a generalized mean-field theory, where we introduce a 
semiclassical distribution function of the center-of-mass 
position of the soliton. Assuming as an initial state an 
isotropic superposition of soliton states, we have argued 
that the measurement process reduces the center-of-mass 
fluctuation and selects a soliton state having a definite 
center of mass. The results obtained by the extended 
mean-field theory and those obtained by the numerical 
simulation using the symmetry-preserving Hamiltonian 
agree excellently in terms of the angular-momentum fluc- 
tuation, the center-of-mass fluctuation, and the conden- 
sate fraction. 

We have thus shown that the critical fluctuations as- 
sociated with the formation of a localized EEC with an 
attractive interaction arise from quantum fluctuations in 
the position of the matter wave, which is also interpreted 
as the depletion of the condensate to the translation 
modes of the bright soliton. 
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APPENDIX A: DERIVATION OF OPERATOR 
RELATION dSp 

In order to prove the relation (|51|l . we disentangle two 
operators 



Oi EE e-i(^^"-^'), 

We define an unnormalized state 
la) =e"*Vac), 



(Al) 
(A2) 

(A3) 



where a is a c number, and define the expectation value 
of Oj {j — 1, 2) with respect to the state \a) as 



(A4) 



The derivatives of Qj with respect to a and a* are given 
by 

^ = {a\d,f^\a) = {a\d,f^d-^d,\a), (A5) 
- {a\TO,\a) = {a\0,Oj'T^O,\a). (A6) 



da* 

By using the relations 

d^^fdi = f coshr-f^sinhr, (A7) 

Oif^O^^ ^ f^coshr + Tsinhr, (A8) 

d^^fOa = f-'q{f + f^), (A9) 

Claf^O^^ = f^-T]{f + f^), (AlO) 

we obtain 

a* coshrQi + sinhr '^^^ , (All) 



da 
dQi 



acoshrQi — sinhr 



da* ' 
da 



(A12) 



da* 

{l-v)a*Q2-V^, (A13) 



da 
dQ2 



(1 - r])aQ2 ~ 1] 



da* 
dQ2 



da* da 
The differential equations for Qj are then given by 

dQi a* + a sinhr 



(A14) 



da cosh r 

dQi a — a* sinhr 
da* cosh r 

dQ2 ^ 1 

da 1 +r] 

dQ2 ^ 1 
9a* 1 + 7] 



Qi, (A15) 

Qi, (A16) 

{a*-i]a)Q2, (A17) 

{a-T]a*)Q2. (A18) 



By integrations of these equations, the expectation values 
Qi and Q2 are found to be 



(A19) 
(A20) 



Q2 OC e 2 l + i? ^^l + rjl"! 2 1 + 7," , 

The operators Oj are thus written as 

0, OC e-W^*'%-'"(i+")^'*e-™^', (A22) 



and we finally obtain 



e 2 



Ivac) = Me 2' 



I vac), (A23) 



e 2 



-S(f-|-ft)2 



|vac)=M2e 2(1+^)*^'' I vac), (A24) 
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FIG. 13: Width of the angular-momentum fluctuations in 
the GP ground state |>I>'''^') (dashed curve) and that in the 



state 6 2 



Crt T )|ygj(,^ (solid curve). 



where the normaUzation constants A/} are determined by 



2?T.! / 1 

1 = (vac|df^di|vac) = AAf ( --tanhr 



(n!)2 \^ 2 



= Wi^coshr, 
1 = (vaclO^iC'slvac) = A/"!^- 
= AA|(l + 77). 



(A25) 



2n! / 1 r/ 



(n!)2 y 2 1 + r; 



(A26) 



APPENDIX B: EVALUATION OF 
ANGULAR-MOMENTUM FLUCTUATION IN 
THE BOGOLIUBOV GROUND STATE 

We justify the approximation, Eq. H54|) with Eq. (|55|) . 
that has been used in order to evaluate the angular- 



momentum fluctuation (A/!^^))^ in the Bogoliubov 
ground state 



(aocj-l- V2ai^^)^|vac).(Bl) 



The angular-momentum distribution of the GP ground 
state l^'e^r, =o) = T^^^o^o + \/2aiS'^)^|vac) is calcu- 
lated to give 



(GP) \|2 



10"^ 10"^ 10" 10' 10^ 10"* lO"* 10^^ 10^ 10^ \{^\'^e,,„,= 



N/2 

E 



1 



"1 



^^{N - C~ 2n)\{C + n)\n\ \ |ao| 

1 £2 

2(AC(«P))2 



4n 



V27r(A/:(GP))2 



(B2) 



where ag and ai are given in Eqs. (jSSl) and and the 
width is given by 



(A£ 



(GP)x2 _ 



2N 



(B3) 



' 2-f |ao|2/|ai|2- 
On the other hand, from the form of Eq. (|52() . the 

angular-momentum-fluctuation operator e^*-^^ •* 
proximately gives the factor 



ap- 



exp 



2(A£(B))2 



(B4) 



to the state with an angular momentum C, where the 
width (A£(^)) is given in Eq. The angular- 

momentum distribution of j^'^^^) is given by the product 
of Eas. llB2|) and (|B4|i . We compare Eqs. (|b5)| and in 
Fig. ^1 where the former is much larger than the latter 
for Ki < 10''. Therefore, the contribution from Eq. ljB2|) 
is negligible, and the angular-momentum distribution of 
can be approximated by Eq. I|B4|I . 
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